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Abstract 

In this article, we consider a Markov process {Xt}t^o, starting from a; € M and solving a stochastic 
differential equation, which is driven by a Brownian motion and an independent pure jump component 
exhibiting state-dependent jump intensity and infinite jump activity. A second order expansion is derived 
for the tail probability ^ x + y] in small time t, for y > 0. As an application of this expansion and 
a suitable change of the underlying probability measure, a second order expansion, near expiration, for 
out-of-the-money European call option prices is obtained when the underlying stock price is modeled as 
the exponential of the jump-diffusion process {Xt}t^o under the risk-neutral probability measure. 

Keywords and phrases: Short-time asymptotics; local jump-diffusion Markov models; stochastic differ¬ 
ential equations with jumps; option pricing. 


1 Introduction 


In this work we consider a Markov process X := {Xt}t^Q with infinitesimal generator of the form 

(1.1) Lf{x) = b{x)f'{x) + ^^-^f{x)+ [ {f{x + -f{x,r))-f{x)-l{\r\^i}-f{x,r)f'{x))u{x,r)dr, 

^ JRo 

where Mq •= ^\{ 0 } b, a, 7 , and ly are deterministic function satisfying appropriate conditions for the 
existence of such a process (see below for further details). Broadly, X can be defined in terms of a stochastic 
differential equation (SDE) of the form: 

dXt = h{Xt)dt + a{Xt)dWt + dJt, 


where W := {Wt}t>Q is a Wiener process and J := {Jt}t>Q is an independent pure-jump process, whose 
jump behavior is dictated by v and 7 as follows: 


IE G [t, f -|- J] : XXg G (a, 6)}] — E [^{s G [t, t -|- J] : G (a, 6)}] 


( 1 . 2 ) 


= E 


r-t+5 I" 

J ^{'rix^.,rMa,b)}i'{X,-,r)drds 


‘Department of Mathematics, Washington University in St. Louis, St. Louis, MO 63130, USA (figueroa@math.wustl.edu). 
Research supported in part by the NSF Grants: DMS-1149692. 

^Department of Mathematics and Applied Mathematics, Virginia Commonwealth University, Richmond, VA 23284, USA 
(yluoSvcu.edu). 


1 



for any t G (0,cx)), 5 > 0, and (a, 6) G ]R\{0}. Intuitively, (1.2) tells us that the jump intensity of the 
process “near” time t depends on its state immediately before t via the function i' in that if is 

large (small), then we expect a higher (lower) intensity of jumps immediately after time t. In the particular 


case of 7 (x,r) = r, (1.2) reduces to 


(1.3) 


E [#{s G [t, t + d]: AXs G (a, 5)}] = E 


rt+S 


viXg -, r)drds 


and u{X^-,r) has the usual interpretation of a stochastic jump intensity as defined in, e.g., [1] and |7j. 
That is, v{x, r) measures the expected number of jumps, per unit time, with size near r when the process 
is at state x. State-dependent jump behavior as described above is an important feature that offers greater 
modeling flexibility to other commonly studied jump processes. For several applications we refer the reader 

to m, m, M, la, 0, m, m, and m- 

The generator 0 covers a wide range of processes. For a Levy processes, h and a are constants, 
7(x,r) = r, and z/(x,r) = h{r), for a Levy density h : E\{0} —)• [0,oo) (i.e., /(x^ A l)h{x)dx < oo). When 
we simply have zz(x,r) = h{r), we recover the class of (local) jump-diffusion models studied in [TO]. In that 
case, X can be constructed as 


(1.4) Xt = x + biX,)ds + a{X,)dWs + Y. 7 iXs-,AZ,)+ Y 7 (^s-,AZ,), 

se(0,t]:|AZ„|>l se(0,t]:0<|AZ„|<l 


c 

where Z is a Levy process with Levy density h and Y^ denotes the compensated Poisson sum of the terms 
therein. The case of z^(x,r) = \{x)p{r) with f p{r)dr = 1 has been studied in |21|. More generally, if 
A(x) := f zz(x, r)dr < oo is locally bounded, we can construct the process X as 

b(X,)ds+ f a{Xs)dWs + Y^iK-^^i)^ 

■^0 n<t 



where 6(x) = b{x) — f l||^|^i| 7 (x, r);z(x, r)dr, 0 < ri < T 2 < ... is a point process on M+ with stochastic 
intensity {A(Xj,-)}s>o and, conditionally on X^- = z, has density ^{z, •)/A(z), independently of any other 
process. 

Unlike the just described processes with finite jump activity (i.e., finitely many jumps during any bounded 
time period), infinite jump activity (UA) processes are important not only from a mathematical but also 
practical point of view. This is especially true in financial applications where several statistical tests, based 
on high-frequency observations, have supported the latter feature. In this work, we study a class of IJA 
processes (i.e., f i'{x,r)dr = oo) that arises as the thinning of a local jump-diffusion process driven by a 
Levy process with stable-like small jump behavior. 


An important obstacle for the application of the class of processes described above arise from the lack 
of tractable transition distributions. In this regard, one stream of the literature has focused on numerical 
methods for the computational simulation of the process (see, e.g., m and m and references therein), which 
in turn can be used to estimate different distributional features via Monte Carlo methods. Another stream 
of the literature has developed approximations for those distributional properties under different asymptotic 
regimes. Approximations in short-time are of particular relevance due to their wide range of applications such 
as statistical estimation and simulation methods. The latter approach was further developed by Figueroa- 
Lopez, et al. m. where a second order expansion, in small time, for the tail probability of the process ( |1.4[ ) 
was developed. Unfortunately, there is almost no work dealing with short-time approximation methods for 
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the distributional properties of state-dependent jump-diffusion. An exception is [23j, where a short-time 
approximation scheme for the transition densities was developed in the special case that ^{x, r) = \{x)p{r) 
with f p{r)dr = 1. 

In this article, we generalize the result in [TO] by developing a second-order expansion, in small time t, 
for the tail probability P {Xt > x + y) (y > 0) of a state-dependent jump-diffusion process X with generator 
(1.1) and initial value x G M. One of the main motivations for considering both the infinitesimal generator 
(1.1) and the tail probabilities is their role in the evaluation of out-of-the-money (OTM) European call option 
prices when the underlying stock price is modeled as the exponential of a state-dependent jump-diffusion X 
under a risk-neutral probability measure: 


(1.5) 


TO = E 



where X has initial value 0 and the log-moneyness k, is such that k > 0. An appealing method for evaluating 

(1.5) is to consider the so-called share measure (cf. [3]), defined as dP^ = e^‘dP. In that case, the 
following neat representation for lit in terms of two tail probabilities holds: 

(1.6) TO = E [(e^‘ - e") lx,>n] = E [e^nx,>.] - e"P [TO >k\=F* [TO > k] - e"P [TO > n] • 

When A is a Levy process, it turns out that the law of X under P^ is again Levy, albeit with different Levy 
triplet, and the same small-time expansion for the tail probability can be used to deal with the two terms in 

(1.6) and obtain an expansion for option price TO (cL jH])- For a general process X with generator (1.1), 
the law of X under P^ is still Markovian with the same generator form (1.1) but replacing u and b with 

(1.7) i/^(x, r) := r), b'^{x) = ^ — j — 1 — e'’'^^’^^l|^l^i 7 (x, r)^ zz(x, r)dr, 

respectively. Note that, in particular, even if the jump dynamics of X were not state-dependent (i.e., 
v{x,r) = v{r)), X would still be state-dependent under P^ (i.e., would depend on r). This motivates us 
to study at once the unifying framework (1.1), which, as explained above, is important in their own right. 

Our main results explicitly quantify the effect of the jump state-dependence in the leading and second 
order terms of the tail probability and OTM call option premium in a small-time setting. Hence, for instance, 
when 7 (x, r) = r (in which case ^{x, r) really measures the intensity of jumps with size near r when the 
process’ state is x), the effect of a positive constant drift b is to change the probability of a large move of 
more than y within time t by 


( 1 . 8 ) 


b(2v{x,y) + J (1-h o(l)), (t^O). 


Similarly, a nonzero constant volatility a will change the probability of a large move of more than y by 

r-oo 


(1.9) 


^ ,2 


. , dv{x,y) 1 

' ^ + 2 


y 


(1 + o(l)), {t 


0 ), 


and the OTM call option premium (1.5) by 

e'’ -|- (9i/(x, r) 


( 1 . 10 ) 




PC 

e^i/(0, k) + 

J k 


dx 


dr + 

x=0 Jk 


(9^zz(x, r) 


dx'^ 


dr 


x=0 


(1 -|- o(l)), {t —> 0), 


which extends a result in na for exponential Levy models. 
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Our approach to obtain the expansion of the tail probabilities of X follows along the same lines of m 
building on a small-large jump decomposition of the process X, similar to that introduced by |18j . An 
essential ingredient of this method is to show that the small-jump component is a diffeomorphism, which 
presents some interesting subtleties in the considered model. More concretely, to obtain the latter property 
under the new state-dependent jump structure, we prove that the state-dependent small jump component 
has the same law as a regular (i.e., state-free) jump-diffusion with sufficiently regular coefficients (see Section 
l^for details). Another novelty in this work is a new, much simpler, proof for estimating the tail probability 
of the so-called one-jump process (the process conditionally on having one big jump at a specified time) 
based on the just mentioned diffeomorphism property, time-reversibility, and a suitable application of an 
iterated Dynkin’s formula (see Lemma 5.1 for details). 


The paper is organized as follows. In Section we formally define the model together with the standing 
assumptions throughout the article. Section [^introduces some needed notations and probabilistic tools. The 
key ingredients to obtain the main result are presented in Sections [^ and The second-order expansion for 
the tail probability is presented in Section while the application thereof to option pricing is developed in 
Section Finally, Section presents a numerical example to illustrate the performance of the expansions. 
To this end, we also develop a simulation method for our state-dependent local jumps-diffusion model based 
on a suitable diffusive approximation of the small jump component of the process. Finally, all the proofs are 
deferred to several Appendix sections. 


General Notation: Throughout, given an Euclidean domain E, C^{E) (resp., C^{E)) represents the 
class of /c*^-times differentiable (resp., bounded and fe*^-times differentiable) functions / : E —>■ M with 
continuous and bounded partial derivatives of order n = 1,... ,A;. In particular, C^{E) C C^{E). Also, di 
and respectively denote the derivative and the n-th order partial derivative operator with respect to the 
i-th variable of a multivariate function. 


2 Setup and assumptions 


In this section, we give a construction of the process of interest and establish the assumptions needed 
throughout. As mentioned in the introduction, we want to consider an infinite-jump activity Markov process 
X with infinitesimal generator (1.1). We use a thinning technique for the construction of the jump structure 
of X based on the jumps of a suitable Levy process Z. To this end, we imposed the following assumption: 


(SI) (i) There exists a Levy density h dominating the jump intensity function i/ : M x Mq (0,oo); i.e., 
v{x,r) ^ h{r)^ for all (x,r) G M x Mq; 

(ii) We also assume that := iy{x,r)/h{r) and h are (^^(M x [€,6]*^) and C^{[e,eY), respectively, 

for any e > 0 and, furthermore, 

liminf inf P(x, r) > 0, limsupsup |r92^'(3^, r)| < oo, limsupsup |5jr'(x,r)| < oo, i = 0,l,2. 

r—> 0 ± a;eK 3;eR r—3;eR 


A necessary and sufficient condition for the first requirement in the previous condition is that h{r) := 
supj, v{x, r) is a Levy density. In that case, we can take h = h] however, in applications it may be more 
convenient to choose another h whose associate Levy process can be simulated more easily. The second 
condition therein imposes some regularity requirements. Note that indeed the derivatives of v do appear in 
the expansions of the tail probability and OTM call option premium (see, e.g., (1.8)-(1.10)), which lead us 
to believe that some smoothness properties on v are needed. 

We are now ready to give the construction for X. Throughout, we consider a filtered probability space 
(n,T', F := {Tijt^OjF) equipped with a standard Brownian motion {Wt}t^o and an independent Poisson 
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random measure p{dt, dr, du) on M_|_ xS := M_|_ xMq (0,1) with mean measure dt h{r)dr du. The compensated 
Poisson measure of p is denoted by p{ds,dr,du) := p{ds,dr,du) — dsh{r)dr du. Then, under the condition 
we have the following construction for the process ;= 

fb{Xi-^)ds+ f aiXi-^)dWs + f [ l\r\>MX^-\r)e{X^Jl\r,u)p{ds,dr,du) 

Jo Jo Jo Je 

+ f f ^\r\^il{^i-\r)0{X^-\r,u)p{ds,dr,du), 

Jo Je 

where 0 : M x Mq ( 0 , 1 ) is a thinning function that takes the form 


( 2 . 1 ) 


( 2 . 2 ) 


9{x,r,u) := 


h(r) f 

Upon the existence and uniqueness of the solution to fc.lL the solution process X^^'^ would be Markovian 


with infinitesimal generator (1.1). The following conditions on 6 , a, and 7 guarantee the well-posedness of 
(2.1) (as proved in Lemma 5.2 below) and other needed features of the process: 

le functions 6 : M —)■ M and cr : M —)• M belong to C'^(M); 


tion 7 (x, r) : M x M —)> M satisfies the following conditions: 

elongs to C^(M x M), and 7 (x, 0 ) = 0 for all x G M; 
all x,r G M, \d 2 'y{x,r)\ > 7 for some constant 7 / > 0 ; 
all X, r G M, |1 + 9 i 7 (x, r)| > ry for some constant ry > 0; 

Some remarks are in order regarding the above conditions: 

iition (S2 ^ is a standard non-degeneracy condition that is also imposed in m- 

iition {S3 ii) implies that, for each e > 0 and i = 0,... ,3, \d\'^{x,r)\ ^ C£\r\ for any |r| ^ e 


and, hence, admits an inverse, denoted by 7 ^(x,r) hereafter. Without loss of generality, throughout 
we assume that r 1 —)• 7 (x,r) is strictly monotone increasing. 


( 2 . 1 ) 

(as j 

(S2) 

(i) 

(ii) 

(S3) 

The 


(i) 


(ii) 


(iii) 

Remark 

1. 

The 

2. 

The 


and 

3. 

The 


4 . Condition (S3 in) is essential for the mapping x —)• x[^'^ to he a diffeomorphism (see Lemma 5.2 below). 


It is worth pointing out that Figueroa-Lopez et al. |10j imposed stronger regularity to the coefficients 
of their SDE. However, we shall see that most results therein are still valid under the milder regularities 
imposed in the present manuscript. 

Remark 2.2. As mentioned in the introduction, in the case that v{x,r) = h{r), for a Levy density h, 
we recover the model (l.f), which was studied in ]10f . Even though it is not evident, under the condition 
(S3 1^, the law of the process (2.1) is actually equivalent to that of a process of the form (I. 4 ) with suitably 
chosen coefficient functions 7 and b (see Remark 4-3 below for more details). However, the resulting 7 
is relatively intractable and does not meet the regularity conditions of m for the second-order expansion 
therein to be applied directly. Furthermore, there are two other important reasons for directly considering the 
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process (2.1). First, the process (2.1) allows the direct modeling and clearer interpretation of the intensity of 


jumps via the parameter u{x,r), which is somehow hidden inside the function 7 in the model (l.j). Second, 
as already mentioned in the introduction, in order to develop the small-time expansion of out-of-the-money 


option prices, one needs to deal with processes having the most general generator (1.1), even if the original 


process X is of the form (l.f). 


Our final assumption is probably the less intuitive. As mentioned in the introduction, there are two key 
ingredients in our approach to obtain the small-time expansion of the tail probability of X. First, we use 
a small-large jump decomposition of the process X. Second, we need that the small-jump component is 
a diffeomorphism. To this end, we use the equivalence of the resulting state-dependent small-jump model 


to a state-free jump diffusion process of the form (1.4), which is possible in light of the above Condition 


as described in the previous Remark |2.2[ However, to conclude that a model of the form ( |1.4[ ) is a 
diffeomorphism, we need some regularity conditions on its coefficients. The main goal of the the following 
relatively mild condition is to establish such conditions (see Proposition |4.1| below). We refer to the Remark 
|4.2| below for further discussion and possible relaxation of this condition. 


(S4) The Levy density h introduced in the Condition (SI) is such that, for some a G (0, 2), g(r) := h{r)\r 
is differentiable in (—eo, 0) U (0, eo), for some eo > 0, and 


Q+l 


liminf g(r) > 0 , 


limsup g{r) < 00 , 

r^0='= 


limsup |rg''(r)| < 00 . 

r^0='= 


3 Some needed notations and preliminary results 

Let Z := {Zt}t^o be a pure-jump Levy process with Levy triplet (0,/i(r) dr, 0), for the truncation function 
The jump measure of the process Z is denoted by q{dt,dr) := j({{t, AZt) G dt x dr : AZt / 0} = 
^{Ti,Ri){dt X dr), where {(Ti, Ri)} are the atoms of the measure q{dt, dr). In that case, the Poisson random 
measure p{dt, dr, du) in ( | 2 . 1 [ ) has the same distribution as a marked point process on {{R, R*)}, with marks 
{Ui}i^i being a random sample from a standard uniform distribution on ( 0 , 1 ). 

In the sequel, the process X^^) is decomposed, in law, into a process with small jumps and an independent 
process of finite jump activity. To formally define the small-jump component of X, we first need to introduce 
a suitable construction for the Poisson random measure p. For any e G (0,1), let 

K{r) ■■= (t)e{r)h{r), he{r) := f)e{r)h{r) := (1 - (j)^{r))h{r), 

where 4>e £ ^“(M) is a “truncation” function such that l^i^e ^ ^ l|r| 7 e/ 2 ) 4>e{'w) is non-decreasing as 

|r(;| increases, and supp(i;()£) = (e/ 2 , 00) U (—00,—e/ 2 ). Next, let Z{e) := {Zt{s)}t^o and Z'{e) := {Zj(e)}i^o 
be independent pure-jump Levy processes, with respective Levy triplets {bz{£), hs{r) dr, 0) and (0, he{r) dr, 0), 
for the truncation function l|r|sgij where bz{e) := r/ie(r)dr. Note that Z{s) is a compound Poisson 
process, and we shall denote its intensity by Ag := f 4i£{r)h{r)dr, and the jump probability density function 
by h^{r) := (j)e{r)h{r)/X^. Let {Nt}t^o, and {Ji}i^i denote, respectively, the jump times, jump 

counting process of Z(£), and an independent identically distributed random sample from the probability 
density function h^. Also, U and J := represent a generic random variable uniformly distributed in (0,1) 
and a generic random variable with the probability density function he{r), respectively. 

The lemma below from [10] will be useful in the sequel. 
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Lemma 3.1. Under the conditions (SI) and (S3) in section^ the following statements hold: 


1. Let 7 (^,r) := z + j(z,r). Then, for each z G M, the mapping r —)■ j(z,r) (resp., r 'y{z,r)) is 

invertible and its inverse ^~^{z,r) (resp., 'y~^{z,r)) belongs to x Mq)- 

2. Both j{z, J) and ^{z, J) admit densities, denoted hyT{r-,z) :=T^{r;z) andT{r-,z) ■.= T^{r\z), respec¬ 
tively, which belong to C^(M x Mq)- Furthermore, they have the representations: 


Feir-,z) = he{j {z,r)) 


d'j 

dr 




-1 


, Teir;z) = hs{j {z,r)) 


d'j 

dr 


(z,7 ^iz,r)) 


-1 


3. The mapping z ^ u := z + 'y{z,r) admits an inverse, denoted hereafter by 'y{u,r), that belongs to 
X 


Now, we are ready to define the “small-jump component” of Let M'{dv,dr) := M'^{dv,dr) de¬ 

note the jump measure of the process Z'{e), and let p'{dv,dr,du) := p({dv,dr,du) (resp. p'{dv,dr,du) := 
p({dv, dr, du) — dvh^{r)drdu) denote the marked point process (resp. compensated marked point process) on 
the atoms of M' with independent uniformly distributed marks on (0,1). For each e G (0,1), we construct 
a process := {Xf (e, defined as the solution of the SDE 


Xf{e,x) = x-\- b^{X(,-{e,x))dv-\- a {X^-{£,x)) dWy 

Jo Jo 

Ns 


(3.1) 


+ E ^ {K- ^ {K- 

i=l 

+ {zil-{^,x),r'^ e (^xl-{e,x),r,u^ p'{dv,dr,du), 

where {lI4}.u^o is a Wiener process independent of p'{dv, dr, du), and 

(3.2) b^{x) ■.= b{x) — / / 'y {x,r) 9{x,r,u)duhs{r)dr = b{x) — / ^ {x,r) v{x,r)4>e[r)dr. 

J\r\!yi Jo Jlr-lsgl 

By comparing their infinitesimal generators, it is not hard to see that the process (|3.1|) has the same dis¬ 


tribution law as the process (2.1) (see Section 2 of [lO] for a more detailed explanation). Next, we let 
X^{e, 0, x) := {Xf (e, 0, be the solution of the SDE: 


(3.3) 


X^(e,0,x) = X-F J bs{Xl-{e,%,x))dv + J cj (^X^_(e,0,x)j dW; 

^ L ^ (^J<:y-{£,^,x),r,u^pfdv,dr,du). 


The law of the process {Xf (e, 0, x)}^^^^^ above can be interpreted as the law of {Xf (e, x)}o<gs<gt condi¬ 
tioning on {Zs(e)}o^s^i not having any jumps. Note that, by conditions (S2) and (S3 i]), the process (3.3) 
is a local martingale with bounded drift whose jumps are bounded by a constant. With equation (9) in |19j 
followed by the proofs of Proposition 3.1 and Lemma 3.2 in [10], we have 


(3.4) 


sup IP[|Xf {rj, 0, x) - x| ^ y] < Ct^ 

0<ri<e,x£R 
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for any y > 0, where > 0 can be made arbitrarily large by taking e > 0 small enough. 

We now proceed to define other related processes. For a collection of times 0 < si < • • • < Sn, let 
{Xf (e, {si,..., be the solution of the SDE: 

{si, ... = X + J be{X^-{s,{si,...,Sn},x))dv + j a(^Xl-{e,{si,...,Sn},x)^dW^ 

+ E ^ [K- ■ • • ’ ^ {K- 

iiSi^s 

+ ^ 0 (^X^-(e,{si ,... ,Sn},x),r,u^ p'{dv,dr,du). 


The law of the process {Xg{e, {si,..., s^}, ®)}o<s<t can be interpreted as the law of {Xf (e, x)}o^s^t condi¬ 
tioning on {Zs(e)}o^<js £4 having n jumps at the times 0 < si < S 2 < • • • < < t. 

For future reference, let us remark that the infinitesimal generator of the small-jump component {Xf (e, 0, x)}t^Q, 
hereafter denoted by L^, can be written as 

Lefiy) ■= T^efiy) +d:j{y) with 


(3.5) 


T^efiy) = beiy)fiy) + ^-^f{y), 

Ief{y)= [ [f{y + liy,r))-f{y)-f{y)'y{y,r)]i^{y,r)^eir)dr, 


where bs is defined in (3.2) and 4>eir) = 1 — ipeix)- Note that, for / G C'^(M), I^f can be written as 


(3.6) ij{y) = j [fiy + 'y{y,rf3)){d2'y{y,rl3)f + f{y + 'y{y,rf3))d^'y{y,r/3) 

-f'{y)dhiy,rP)] {1 - P)dPi'{y,r)(j)e{r)r^dr, 

which is finite and bounded due to the conditions ( STp ) and (S3 I]). 

The following first and second order Dynkin’s formula for the process {Xf (e, 0, x)}i^o will be needed in 
the sequel: 

(3.7) E [/(Xf(e,0,x))l = /(x) + t E [L,/(X^,(e, 0, x))! da, V/ G ^^(M), 

(3.8) E[/(Xf(e,0,x))l =/(x)+tL,/(x)+t2 /" {1 - a)E\{L,f f{xUe,(l},x))\da, V/G Q^(M). 

Furthermore, for / G C'^(E) (resp., / G (7^(M)), sup^. |L£/(x)| < oo (resp., sup^, |Lg/(x)| < oo) and, thus, 
the reminders in (3.7) (resp., (3.8)) is 0{t) (resp., 0{t^)) uniformly on x and t. The proofs of (3.7)-(3.8) 
follows from Ito’s formula and goes along the same lines as in the proof of Lemma 3.3 in |10j . 

Remark 3.2. By the condition (S3 ^ in Section^ the mapping r —)■ 7 (x,r) is a bijection in M. Then, the 
generator (3.5) can be rewritten as 

(3.9) LJ{y) = be{y)f{y)+ ^-^fiy) + j {f{y + z)-f{y)-zf'{y))K,{y,dz), 
where the kernel K^{y, dz) can he written in the following two equivalent forms for any A G 13{R.), 

(3.10) Ke{y,A)= [ lA{ 7 iy,r))i^{y,r)^eir)dr = [ n {y,-f-^{y, z)) {l~^iy, z)) dz-f~^{y, z)dz. 

dR Ja 
























4 A weak solution process 


The main purpose of this section is to introduce an approach to overcome the difficulties posed by the 


discontinuous jump component 'y{x,r)8{x,r,u) in (3.1) and (3.3). To this end, we first show that the state- 
dependent local jump diffusions (3.1) and (|3.3|) are equivalent to a state-free local jump-diffusion of the form 


(1.4) and prove some needed regularity on its coefficients. 


Note that the process {Xf{e, 0, x)}t^o is a semimartingale (see, e.g., [T5l III.2.18]) and, furthermore, com¬ 


paring the generator (3.9) to that in m IX.4. 6 ], it is a homogeneous diffusion process with jumps as defined 
in |15L IX.4.1]. Then, we can determine the semimartingale characteristics {B*,C*, 1 }*) of {Xf (e, 0, 
relative to the identity truncation function, as 


/'t ci 

( 4 . 1 ) B^= (e,0,ds, = / cr^ ^Xf (e,0,x)^ ds, 'd*(dt,dr) = Xg ^X®_(e,0,x),dr^ 


dt. 


By Definition III.2.24 and Theorem III.2.26 in [15], a semimartingale with characteristics specified by (4.1) 
is a weak solution of the SDE 

(4.2) Yt{£,^,x)=x+ f be{Ys{£,<I},x))ds + f a{Ys{e,%,x))dW* + f f 6{Ys-{£,<I},x),r)]I*{ds,dr), 

Jo Jo Jo J 


where, under the solution measure P* in the canonical space (D*,X*,P*), {W*}s^o is a one-dimensional 
Wiener process, fi*{ds,dr) is an independent Poisson random measure on x Mq with intensity measure 
dsF{dr) and corresponding compensated measure Jl*{ds,dr) := fj.*{ds,dr) — dsF(dr). Here, F is a positive 
(T-finite measure on Mq to be chosen below, while d : M x Mq 1^ is a Borel function implicitly determined 
by K, 5, and F via the equation 


(4.3) 


Ks{x,A) = J lA(d(x,r))F(dr), e H(Mo), 


where we recalled from (3.10) that K^(y,A) := J^lA{'y{y,r))i'{y,r)(ps{r)dr. In what follows, we take 

F{dr) = ^e{r)h{r)dr, 


for e < eo, with eo given as in Condition (S4). 


In order to identify the function d(x, r) corresponding to the above measure F, we introduce the following 
functions : Mq 1^0 and ■0 : M x Mq 1 ^ 0 - 


(4.4) = < 


poo 

/ j)s{r)h{r)dr, w > 0, 

J W 

(f)e{r)h{r)dr, w < 0, 


^p{x, w) = < 


f°° 

— / (j)s{r)i'{x,r)dr, w > 0, 

J W 
w 

^s{r)i'{x,r)dr, w < 0. 


Note that the restrictions of the two mappings w 1 —and w 1 —)■ 'tl>{x,w) on := (—£,0) U (0, e) are 
strictly increasing and one-to-one, with range Mq- Thus, they admit “local” inverses defined on Mq with 
range D^, hereafter, denoted by 'iIj~^{x,w) = ip~^{x,w) and 'ip~^{x,w) = '4)~^{x,w), respectively. 

) and the fact that the support of the mean measure ds(l)^{r)h{r)dr of y* 
lies in (—e,e), it is clear that the values of 6{x,r) for jrj > s are superfluous and we only need to define 
d(x.r) for r G D^. Let d : M x —>■ M be dehned as 

(4.5) d(x, w) := 7 (x, '0“^(x, 'tp{w))) . 


From the form of the SDE (4.2 
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In order to show that the above function satisfies (4.3), observe that, for each x G M and w G L>e, the 
mapping w —)> 6{x,w) is strictly monotone and, thus, its inverse, hereafter denoted by 6~^{x,w), exists and 
satisfies 

(4.6) 6~^{x,w) = ijj~^ ['ip{x,'-f~^{x,w))) . 

Then, ^p{6~^{x,w)) = 'tp{x,'y~^{x,w)) and, from the definitions of V' and we have the identity 


JS ^{x,w) 
rS~^ (x,w) 


cj)s{z)h{z)dz, w > 0, 


= 


(f>e{z)h{z)dz, w < 0, 


- / 4>£{z)iy{x, z)dz, w > 0, 

p'y~'^(x^w) 

/ (j)i;{z)iy{x, z)dz, w < 0. 

' —oo 


Upon differentiation with respect to re, we get 

(4.7) (t)e{5~^{x,w))h{5~^{x,w))d25~^{x,w) = ^e{l~^{x,w))u{x,-f~^{x,w))d2'y~^{x,w), 


which implies (4.3) with the chosen measure F{dr) = (f>^[r)h{r)dr. 

We now proceed to show some needed regularity properties of the function 5 : M x 


i, with which the 


almost sure existence of a stochastic flow of diffeomorphisms associated to the SDE (4.2) below is guaranteed 
by virtue of results in [3] (see Lemma 5.2 below). The proof of Proposition 4.1 is deferred to the Appendix 

EH 


Proposition 4.1. Under the conditions (SI), ^83^ , and (S4) in Section^ the function 6 defined in (4-5) 
can he continuously extended on M x (—e,e) with (5(x,0) := 0, for any e > 0. Furthermore, for e > 0 small 
enough, 

(4.8) (i) sup |92 (^(x, r(;)| < oo, (ii) |cl^(5(x,t(;)| < (hi) |1 + 5i(5(x,t(;)| > ??, 

for any 0 ^ i ^ 2, w £ D^, and some constant r],k > 0, independent of x and w. 

Remark 4.2. The main purpose of the assumption (S4) in Section^ is to simplify the verification of the 
regularity of 5 as stated in Proposition 4-1. However, it is important to note that all what follows, including 
the main Theorem 6.1. hold true if u and h are such that the function F{x,w) = fi~^{x,ip{w)) satisfies the 
following conditions: 


(4.9) 


sup |cI^F(x,ro)! < oo, sup \dl,F{x,w)\<oo, i = 0,1,2. 


In particular, if iy{x,r) = h{r), for an arbitrary Levy density h that is smooth enough outside any neighbor¬ 
hood of the origin, then trivially holds true since F{x,w) = w. In that case, we recover the results in 

m- 

Remark 4.3. Using similar arguments to those at the beginning of this section, it is not hard to check that, 
under the condition (SS'^, the model (2.1) is equivalent in law to the model (1.4) with 7 replaced by an 
appropriate function T. Concretely, we need to take 


T{x,w) ='j (x,fi ^ (x,'fi{w 
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where 


w > 0, 


ip{w) := < 


POO 

/ h{r)dr, 

J W 

/ w 

h{r)dr, w < 0, 

-OO 


tpix, w) := < 


v{x, r)dr, w > 0, 

V 

z^(x, r)dr, w < 0. 


However, the regularity o/T is harder to study than that of the function 5 in (4-5). Indeed, for instance, the 
first-order partial derivatives of the function F{x,w) := 'ilj~^{x, if{w)) are given by 


d2F{x,w) = 


h{F{x, w)) 


diF{x, w) = — 


{di'ip){x,F{x,w)) 


v{x,F{x,w)) v{x,F{x,w)) 

and the behavior of h{w) and v{x, w) as w ^ ±oo is now relevant as well. 


5 Tail estimate for the one-big-jump process 


In this section, we give an expansion for the tail probability of the process defined in (3.1) conditioned 
on Ns having only one jump. The following lemma (proof in Appendix ] A. 2| ) is a counterpart of Lemma A.l 
in m, even though the proof developed in the present article is new and much simpler. Below, we recall 
the notation v{x) = (T^(x)/ 2 and set 

ne{x,r) := u{x,r)(j)e{r), iy^{x,r) := i/(x,r)(l - feir))- 
Lemma 5.1. With the notation introduced in Section^ let 


(5.1) 


H{t, z, q) := E 


v{z,J) 


Af(e,0,u) ^ q 


v=z+-^{z,J) 


V HJ) 

Then, under the conditions (S1)-(S4) in Section^ 

(5.2) H{t, z, q) = Ho{z; q) + tHi{z; q) + t^izf\z, q), 

where. 


Ho{z;q) 

D{z-q) 

I{z-,q) 


— / Ve{z,r)dr, Hi{z;q) ■.= D{z]q) + I{z-,q) 

^([^£( 9 ) iz,T^{z,q)) d2T^{z,q) 

- v{q)[d 2 i'£ {z,j~^{z,q)) {d 2 T^{z,q))‘^ + {z,j~\z,q)) dlT^{z,q)\^, 


1 


I ^ 

/ W {z,j~^{z,r])) d2j~\z,ri)ue{v,r)dri 
.J'yiq,r) 

- Ve (z, 7 “^(z,g)) d 2 T^{z,q)-i{q,r)h'e{q,r) dr, 


and, for e > 0 small enough. 


(5.3) 


lim sup sup 

t—^-O z,q 


(z, q) 


< OO, 


sup \Hi{z-, g)| < 00 . 

z,q 
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The proof of Lemma |5.1t which can be found in Appendix A.2 builds on the following key lemma. The 
proof can be found in Appendix |A.3 


Lemma 5.2. Under the conditions {S1)-{S4) in Section^ the SDE [4.2) admits a unique solution, which in 
turn implies the existence of a unique weak solution of the SDE (2.1). Moreover, for any t > 0, the mapping 
X I— )■ yj(e, 0,x) is a diffeomorphism on M. 


6 The second order expansion for the tail probability 


In this section, we obtain a second order expansion, in a short time t, for the tail probability P X): ' ^ x + y 

for any y > 0. The idea is to exploit the equivalence in law of and the process X^ defined in (3.1), and 
to condition on the number of jumps of Z{£). Specifically, we have 


( 6 . 1 ) 


X^""^ ^ x + y 


Xf(£,x) ^ x + y = ^ P Xf (e, x) ^ x + 


Nf = n 




n=0 


n\ 


In order to present the expansion, let us hrst recall the notation u^{x,r) := i'{x,r)(j)c{r) and Usix,r) := 
v{x,r)(()e{r) as well as the functions 7 , 7 ”^, 7 ~^, and 7 , introduced in Lemma 3.1 The following theorem 
(proof in Appendix |A.4[) states the main result of this article. 


Theorem 6.1. Under the conditions (S1)-(S4) in Section^ the following asymptotic expansion holds, for 
any y > 0, 


( 6 . 2 ) 


^ x + y 


= tPi{x,y) + —P 2 {x,y) + o(t^), 


as t ^ 0, where the coefficients Pi{x,y) and P 2 {x,y) admit the following representations, for e > 0 small- 
enough: 

Pi{x,y):= i^eix,r)dr, P 2 {x,y) := Ve{x,y) + Je{x,y), 


with 


T^e{x,y) = beix) 

(T^(x) 


+ 


roo 

-ne {x,'y~^ix,y)) [di'y~^{x,y) - d 2 'y~^ix,y)] + / di^eix, r)dr 

J'y-^(x,y) 

2 

d 2 T^e{x,l~^{x,y)) {di'^~^{x,y) - d 2 'y~^{x,y)) + dfi's{x,r)dr 


-Ve{x,-/ ^{x,y)) [ 5^7 ^{x,y) -28182^ ^(x,y) + a |7 
+ {be{x + y) - a{x + y)a'{x + y)'^Us (x,7"^(x,y)) 82')~^{x,y) 


a‘^{x + y) 


{d 2 Ve{x,-i ^{x,y)) [027 ^(a;,y)] (x ,7 ^(x,y)) 9|7 ^(x,y)), 
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Je{x,y) = 


/ i^eix+ ^{x,r),ri)dri - i^s{x,ri)dri 

+ 'y{x,r)(i^e{x,'y~^{x,y)) [di'y~^{x,y) - d 2 'y~^ix,y)] - / diUeix,ri)dri 


I'eix, r)dr 


+ 


fx+y 


i^e{x,^ ^{x,ri))d2j ^{x,ri)iyeiri,r)dri 


J L J‘y{x-\-y,r) 

- {x,-f~^{x,y)) d2'y~^{x,y)-i{x + y,r)i?e{x + y,r) dr 

/ poo 

Ve{x,ri) / V^{x,r2)dr2dri 

J 7“i (a:+7(x,ri),y—7(a:;,ri)) 

poo p poo p 

- / i/£(x,ri) / Us{'y{x,ri),r2)dr2dri - / Ve{x,r)dr / Ve{x,r)dr. 

J 'y~^{x,v') J J ''i~^{x.v') J 


Remark 6.2. 


1. The expansion (6.2) indeed does not depend on h, which is expected since the infinitesimal generator 
(1.1) of X only depends on ( 6 , a, 7 , z^). Obviously, the coefficients Pi and P 2 are independent of e, 
even though, the given representations involve e. In particular, note that, since 7 (x, 0 ) = 0, for s > 0 
small-enough, Pi{x,y) = J^^.^^^.^^y,^iy{x,r)dr. 

2. The drift and diffusion are absent in the leading term, which can be interpreted by saying that a possible 
“big” move of the proeess in a short time happens mostly as a result of a “large” jump. This phenomenon 
also appears in the expansion of m- In the case that u{x, r) is interpreted as the probability density 
of a mark J of the underlying marked point process driving X, then Pi{x,y) = P (y{x, J)fiy]. 

3. It is tedious but not hard to verify that the above expansion reduces to the expansion of m when 
v{x,r) = h{r), in which case the jump intensity does not depend on the state of the process {Xt}t^o. 

4- When 7 (x,r) = r and e G (0,1) is small enough, the expansion reduces to 


t / v{x,r)dr 


+ 7 I be{x) 


v{x,y)+ / div{x,r)dr 
Jy 


a^x) 


1*00 

02 ^ (^,2/)+ / dfu{x^r)dr 

Jy 

+ {be{x + y) - a{x + y)a'{x + y)'^v {x,y) - ^ ^ 


2 

/ roo poo / poo 

/ z^(x + r, ri)dri — / z^(x, ri)dri — r ( z/(x, y) + / div{x,ri)dri 

J v—r Jv \ Jv 


-d2v{x,y) 


\_Jy-r 


+ f [ ne{x,ri)nfix + ri,r)dri-ns{x,y)rnfix + y,r) 
J J V—r 


LJy-r 


dr 


z^e(x, r2)dr2dri 


poo 

(x,ri) / 

ly-ri 

i^(x,ri) J ^^{xri,r 2 )dr 2 dri — J v{x,r)dr J z/£(x, r)dr| + o(t^). 


t'£(x, r)dr 
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In particular, supposed that b and a are constants. Then, recalling that = b — J^^^^.^rn{x,r)4>e{r)dr, 
the effect of a positive “drift” b is to increase the probability of a “large” move of more than y by 

ff ( d \ 

—b\^2u{x,y) + J —u{x,r)drj{l + o{l)). 

Note that the second term inside the parentheses above is missing in the absence of a state-dependent 
jump intensity. Similarly, a nonzero constant volatility a will change the probability of a “large” move 
of more than y by 


2 


a 



v{x,y) + 


1 

2 



dx'^ 


u{x, 



(1 + o(l)), 


in short-time. Again, the seeond term inside the parenthesis is the effect of a state-dependent intensity. 
For a general b, it is also intuitive that the net effect depends on the average drift, {b{x) + b{x -\- y))/2, 
at the initial and final points x and x + y. The net effect of a general function a on the probability of 
a large positive move of more than y depends on the functional 


Hx,y) 


A 

dy 



+ cr^(x + 


+ 





as ^A{x,y) (1 + o(l)). 


7 The small-time second-order expansion for OTM call prices 


In this section, we derive a second-order expansion, in short-time, for the price of an out-of-the-money 
(OTM) European call option, with maturity t and strike K, written on a nondividend paying stock, whose 
risk-neutral price process is modeled by 

v(o) 

St = Soe^t , t ^ 0, 


where the process {xj:^'^}t^Q is given by (2.1) with the initial condition x 
this section, we omit the superscript in XOb 


0. For simplicity, in the rest of 


As explained in the introduction, we shall consider the share measure associated with the stock (namely, 
the martingale measure obtained by taking the stock as numeraire) to evaluate the premium of the OTM 
option. Concretely, let k := log{K/So) be the so-called log-moneymess of this call option and as customary 
suppose that the risk-free rate r is 0. Then, the price of this option can be written as 


(7.1) E[(Si - K)+] = E[(Si - K)1s,;,k] = [Xt ^ k] - Soe’^F [Xt ^ k ], 


where is a probability measure, locally equivalent to E, defined as E'^[i?] := E[e'’^‘l 5 ] for any T)- 
measurable set B. Hereafter, denotes the corresponding expectation. For E^ to be well-defined, 
{St/Sojt^Q must be a F-martingale, for which we impose the following drift restriction; 


(7.2) 


b{x) -I- 


a\x) 


_(_ / - 1 - l|r|<gi 7 (x, r)) i^(a:,r)fir = 0. 


The integral in (7.2) will be well defined under the conditions 
condition: 


and (SI) in Section as well as the 



g7(a:,r’) 


i>{x, r)dr < oo. 
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That {Xt}t^Q is a F-martingale under (7.2) is a consequence of Ito’s formula. 

The second-order expansion for the probability appearing on the second term of (7.1) was treated in 
Section]^ Next, we seek for a second order expansion in t for the tail probability ^ k]. To this end, 

we impose the following condition: 


(S5) The function g introduced in Assumption (S4) is such that e^'^g{r)dr < oo, where c is defined by 
c := sup | 927 (x, r)| < oo. 

x^r 

Our hrst task is to determine the infinitesimal generator of the process under F^. To that end, 

we compute the expectation 

E# [q{Xt)]=¥.[e^^q{Xt)] 

for an arbitrary function q G C'^(M). Let f{x) = e^q{x). Then, applying Ito’s formula (see, e.g., [T] Theorem 
4.4.7), / {Xt) = Mt + At, where 

Mt = q{0) + f [q (A,) + g' (A,)] a (A,) dW, 


'lO,t)xE 


’^)q{Xs- +'y {Xs-,r)) - q{Xs-) lg(^^^_^^^^^^^p{ds,dr,du), 

At = f [q (A,) + q' (A,)] b (A,) ds+l f [q (A,) + 2q' (A,) + q" (A,)] ^7^ (A,) ds 

Jo ^ Jo 

eA^‘’Aq (A^ 7 (A*, r)) - q (A^) - l|r-|^i 7 r) [q i^s) + q (A^)] ^ u (A^, r) drds. 


+ / e 
Jo 


Thus, 


where 


¥*[q{Xt)]=q{Q)+ / E# L*q{X,) 


ds, 


L*q{x) = b{x) [q{x) + g'(x)] -h [g(x) -b 2q{x) + q"{x)] 


-b 


eA^’Aq{x + 7 (x, r)) - q{x) - l|r|^i 7 (x, r)[q{x) + q\x)] v{x, r)dr. 


Comparing the above formula to the Dynkin’s formula (3.7), we can then identify as the infinitesimal 
generator of {Xt}t^o under F^. Using the martingale condition (7.2), we can further write as 

(7.3) L*q{x) = b*{x)q{x)-\ -^Q'"(x) + / [q'(x + 7 (x, r)) - g(x) - l 7 |^i 7 (x, r)g'(x)] i^^(x, r)dr, 

where 

(7.4) .7^.0 := b*(.) := MU + U(.) + / _ l) .)*. 

Note that under the Conditions (S3 I]) and (SI) in Section]^ the integral appearing in (7.4) is well-defined 
and, furthermore, it is not hard to see that belongs to C? and can be written as 


(7.5) 


b*{x) = ^ - I (eA-x) _ 1 _ u{x,r)dr, 
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in light of (7.2) 


In order to conclude the second-order expansion for the call option price, it remains to show that the 
measure satishes the Conditions (SI) and (S4). This is obtained in the following result whose proof is 
deferred to Appendix |A.5 


Corollary 7.1. Under the Conditions (S1)-(S5), for any y > 0, we have 


[Xt ^y]= tP*{0, y) + (0, y) + o{t^), 


as t —?• 0, where, Pf" and P^ are given as Theorem 6.T but replacing v and b by and 6^, respectively, as 
defined in ( 7.4)- 


Finally, using the pricing formula for the European call option introduced at the beginning of this section, 
the price of an OTM European call option has the following expansion in a short maturity t, for A; > 0: 


E (5* - Soe^^ = ^ kj - [Xt kj 


(7.6) 


= tSo[Pf(0,k)-e^Pi(0,k)]+jSo[P2*(0,k)-e^P2(0,k)]+o(t^). 


The formula (7.6) extends the hrst order expansion for the option price given in |10j for a non state-dependent 


jump intensity process (i.e., when v{x,r) = h{r)). 

Remark 7.2. 


1. It is not surprising that the leading term is only determined by the jump component according to the 
formula 

(7.7) tSo ( [ e'^^^'‘^^i^{0,r)dr — e^ f n{0,r)dr \ = tSo [ (e"^^^’^'^ — e^'] ufi),r)dr. 

{r:"/{0,r)^k} J {r:'y{0,r)^k} J J ^ 

In particular, if C{t, K) := E [St — AT) , represents the premium of a call with expiration t and strike 


K and 'y{x,r) = r, ( 7.6)-( 7.7) suggest that, for k := log{K/So) > 0, 
d^C{t,K) 


dK^ 


te 


.(0 


thus, the curvature of the call premium K ^ C{t, K) is strongly determined by jump intensity n{0, k). 


2. Plugging the expansion given in Remark 6.2-4 'I'^to ( f7.6| ) and using (7.2) and (7.4), we note that, when 
7 (x,r) = r, the effect of a nonzero constant volatility cr(x} = a in the price of an OTM call option is 
of order 


So 


t^a^ r 


eV 


/■oo r gfc rc 

(0,A:)-F / —-— divfi),r)dr + 

Jk ^ Jk 


e — e 


dfiy(0, r)dr 


(1 -|- o(l)), t —> 0, 


which extends a result in /71]/ for exponential Levy models. We again observe an additional contribution 
due to the state dependent feature of the model. 
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8 A numerical example of tail probability estimates 


In this section, we show a comparison of numerical estimates for the tail probability ^ y] by the first 
and the second order expansion in a small t, and a Monte Carlo estimate based on the simulation of {Xt}t^o 
with a jump augmented Euler-Maruyama scheme (cf. |21)1. combined with a diffusion approximation of the 
small-jump component of X. For the numerical results, we use the following parameters 


( 8 . 1 ) 


b{x) = sin X, a{x) ~ 2 4 ~ 

/ 1 3 

h{r) = {a = 1 . 01 ), r) = ( —tan~^{x) + - 

V ZTT 4 


I — 1—a 


=: c{x)h{r), 


which satisfy the conditions (S1)-(S4) in Sectionj^ For simplicity, we set = l{r:|r|>£}) when computing 
the coefficients Pi{x,y) and in the expansion (6.2). This is valid since Pi{x,y) and P 2 {x,y) don’t 


depend on and, thus, one can consider a sequence of smooth truncation functions, that converges to 

0£,(r) l{r:|r|>£}? US Tl ^ CX). 

8.1 Second order approximation for the tail probability 


With the parameters (8.1) above and fixing x = 0, we can compute explicitly the expansion stated in Remark 


6.2-4, in which 


7 ( 0 ,r) = 7 ( 0 ,r) = r, 7 ( 0,?/)=7 { 0 ,y)=y, u^{ 0 ,'y { 0 ,y)) =-he{y), 

di 7 ( 0 ,y) = di 7 “^( 0 ,y) = djd^ 7 ( 0 ,y) = dj^ 7 "^( 0 ,y) = 0 {i + j = 2 ), 




d27(0,y) = d27 (0,y) = d27 (0,7/) = 1, Ml{ 0 ,ri),r 2 ) = c{ri)heir 2 ), d2i^s{0n {0,y)) =-h^{y). 


Then, 


-Pi (0,7/) = 


hs{r)dr, 


'i^e{o,y) = -^Kiy) + ^h^{y)[be{y) - (^{yWiy)] - ^o-^(y)^e(y)> 


X(0,7/) = - 


1*00 q poo / q 1 /*00 

z{r) / heiri)dri - - j h^{ri)dri - r ( -K^y) + — 

J v—r 


27r 


poo 

/ he{ri)dri 
Jv 


^3 


c(ri)/i£(ri)dri - c{y)he{y)r 


L Jy-r 


h^{r)dr 


9 f f°° 3 f°° 9 f°° 

^ 16 J J hs{ri)dridr --X^ J he{r)c{r)dr -—X^ J he{r)dr. 


8.2 A Monte Carlo estimate of the tail probability 


he{r)dr 


In the sequel, we present a numerical method to simulate the process X := X^^'^ defined in SDE (2.1). This 


is based on a diffusion type approximation of the “small-jump” component of the process together with a 
jump augmented Euler-Maruyama scheme (cf. [ 2 T]) for the “big-jump” component. To introduce the main 
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idea, let us start by writing the infinitesimal generator L of X, defined in (1.1), as 


Lf{x) = b{x)f{x) + 


( 8 . 2 ) 


+ 

+ 


f{x + 7 (x, r)) - fix) - 7 (x, r)/'(x)l||.^(,,^^)|^i} uix, 
fix + 7 (x, r)) - fix) - fix)-fix, r) i^(x, r)l{\^(^^^r)\^^}dr, 


for any / G C|(M) and e G (0,1), where 5(x) := 6(x) — f 7 (x, r)(ir. Note that, 

as e —7- 0, the term in (8.2) can be approximated by 


Therefore, for e small, the generator L is “close” to an operator on C|(M), dehned as 

(8.3) LJix) :=bix)f'ix) + ^^^^f"ix) + J fix + z)-fix)-zf'ix)l{\^\^iy K^ix^dz), 

where 

alix) := a‘^ix) + alix), Ksix,A):= j lAi'yix,r))i^ix,r)l{i^^^^r)\>e}dr. 

The operator is the infinitesimal generator of a Markov process {Xtie, x)}t^o of finite jump activity. 
Indeed, as stated in Remark |7.2[ condition ( S^ ) implies that there exists an eo > 0,^depending only on e 
and not on x, such that {r : |7(x, r)| > e} C {r : |r| > eo}) for all x G M. Therefore, {Xf}t^o can be defined 
as the solution of the SDE: 


(8.4) 


Xf = x 


t t 

+ f kiX^^.)dv+ f + 

Jo Jo 


where beix) := bix) - / 7 (x,r)[l||r|,gi} - l{|^( 3 ,_r)|^£}]z/(x,r)dr, %ix,r) := 7 (x,r)l||.^(a._^)|>£}, is 

a Wiener process, is an independent Poisson process with jump intensity = / hir)l^^j.\>s}dr 

and jump arrival times are independent identically distributed with the probability density 

function hdr) = /i(r)l||j,|>£}/Ae, and is a random sample from a standard uniform distribution. 

The following result, whose proof is deferred to Appendix A.6, justifies that can be asymptotically 

approximated in law by {X^}t^Q, as e —)■ 0. 


Lemma 8.1. 


Under the conditions (S2)-(S4) in Section ^ {Xf}t^Q ^ {J^t} 


!t>o- 


8.3 Implementation and numerical results 


By Lemma 8.1, we can approximate {Xt}t^o by simulating the process defin ed in (8.4). Since 


{m is a Poisson process with intensity A^, the inter-arrival jump times of ( |8.4[ ) have independent 
identical exponential distribution with mean 1/Ag. Then, we generate the jump times {ri,r 2 , • • • ,Tm} by 
setting 

1 

To = 0, Tk = Tk-l + T-efc, 
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where m := max{A; ■ Tk ^ t} and {ek}k>o is a random sample from the exponential distribution with 
parameter 1. Also, we use the inverse transformation sampling method to obtain a sample of the jumps 
{Jj}™ Finally, we construct over the jump-augmented time steps 

{<UEr u {t.}™ „ = A n € N+) , 

The algorithm is similar to the one in [21] but with an additional thinning condition at each jump time r,. 

Figures |8. If Figure |8.3| show the tail probability ¥[Xt ^ y] estimated by the first and the second order 
expansions and the Monte Carlo approximation, for e = 0.1,0.01,0.001. We see that the second order 
expansion is indeed a better estimate than the first order expansion. We also observe that as e —)• 0, the 
second order expansion approaches to the Monte Carlo approximation and, furthermore, that the latter 
exhibit little variation for different e, which in some sense justifies the diffusion approximation of the small- 
jump component of X. 


Tail probability estimates (s=0.1, t=10’'^) 



Figure 8.1: Expansions and Monte Carlo approximation of P[W ^ y\ with e = 0.1. 
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Tail probability estimates (s=0.01, t=10"^) 



Figure 8.2: Expansions and Monte Carlo approximation of P[Xt ^ y] with e = 0.01. 


>• 10 ^ Tail probability estimates (s=0.001, t=10"^) 



Figure 8.3: Expansions and Monte Carlo approximation of ^[Xt ^ y] with e = 0.001. 


A Proofs 

A.l Proposition |4.l] 

Proof: The first assertion holds true since lim^_^o± = Too and lim^_^±oo '4’~^{x,w) = 0 and, thus, 

lim S{x, w) = j I X, lim (x, ip{w)) ) = 7 (x, 0) = 0. 
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This, shows that the function 5 can continuously extended on M x (—e, e) by defining (5(a:,0) := 0. Now, we 
proceed to show the first two assertions in ( |4.8[ ) for any w G (0, e). We can similarly show the case w G (—e, 0). 
Throughout, we use the notation F{x, w) = {x, ^^{w)) and recall that 5{x, w) = 7 (x, F{x, w)). For future 

reference, let us also note that, for any w G 

(AT) |F(x,u;)| < Ire] and (j)e{w) ^ ^e{F{x,w)), 


because 

roo roo pw pw 

/ (^e(r)h(r)dr ^ / ^^[r)v{x,r)dr {w > 0), / ^^{r)h{r)dr ^ / cj)^{r)v{x,r)dr (u) < 0), 

J w J w J —oo J —oo 

and the function ^e{w) is chosen to be decreasing as |i(;| increases. 

We start to prove (4.8-i) for w G (0, e). To this end, note that dwd{x,w) = {d 2 'y){x, F{x,w))dwF{x,w) and, 
by Condition ( S^ ), it suffices to show that \dujF{x,w)\ is bounded. By the definitions in Eq. ( |4.4| ) and 
recalling that v{x,w) = v{x,w)h{w) and h{w) = g{w)\w\~°‘~^, 


(A.2) 


dwF{x,w) = 


ijj'lw) 


9{w) 


{d 2 'ip){x,F{x,w)) u{x,F{x,w))g{F{x,w)) \ w 


F{x, w] 


Q+l 


Mw) 


(l)e{F{x,w)y 


for any w G (0,e). Note that, due to ( A.l| ) and Conditions (SI ^ and (S4), for e > 0 small enough, 
\rF.xi>{x,F{x,w)) > g and infa,g'(F(x,rc)) > g, for any 0 < w < e. Using the latter and again (A.l), we 
conclude that 

sup \dwF{x,w)\ < oo, 

0<w<e,x^M. 

which, as discussed above, implies ( |4.8| -i). 

We proceed to show (4.8-ii), for i = 0,1, 2. The case of i = 0 follows from (4.84) by the mean value theorem 
and the fact that 6{x, 0) = 0. For i = 1, note that dxS(x, w) = {di'y){x, F{x, w)) + {d 2 g){x, F{x, w))dxF{x, w). 
By Condition (S3 i]) and (A.l), 


|(ai7)(x,F(a:,u;))| < 


sup \d 2 di'y{z,r)\ 


\F{x, re)! < K\w\, 


for some constant K. Similarly, since sup^^,, |cl 27 (^, ?’)| < oo, it suffices to show that \dxF{x, w)/w\ is bounded. 
To this end, note that 


{diy){x,F{x,w)) _ fp^^^^^diu{x,r)g{r)r ^ Vg(r)dr 

{d2'4)){x, F{x, w)) P(x, F{x, w))g{F{x, w))^e{F{x, w)) 


Since dii'(x,r) and g{r) are bounded for r in a small neighborhood of the origin and, again, (j)e{w) is 
decreasing in |r(;|, for some constants K, K' < oo and any 0 < w < s, 


dxF{x,w) 


w 


< K 


F{x,w)°‘+^ " ^(j)e{r)dr 

wi>{x, F{x, w))g{F{x, 'w))4>s{F{x, w)) ~ i>{x, F{x, w))g{F{x, w))^siF{x, w^w ’ 


< K' 


F{x, wy~^^(f)s{F{x, w))F{x, w) 


and one can use the same arguments as those in (A.2) to show that \dxF{x,w)/w\ is bounded on M x (0,e) 
and, hence, to conclude the validity of (4.84) for i = 1. For future reference, note that the previous arguments 
also show that \dxF{x,w)/F{x,w)\ is bounded on M x (0,e). 
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For i = 2, note di5{x,w) can be decomposed as 


{dfj){x,r) + 2{d2di'y){x, r)dxF{x,w) + (9|7)(x,r) {dxF{x,w)Y + {d 2 'y)ix,r)dlF{x,w) 


r=F{x,w) 


The first three terms are trivially bounded by k\w\, for some constant k, by Condition (S3 (A.l), and the 

boundedness of \dxF{x,w)/'w\. It remains to show that \d‘^F{x,w)/'w\ is bounded. From straightforward 
differentiation, 


dxF{x,w) = - 

+ 


{dlY)ix, F{x, w)) _ 2 {d 2 diY){x, F{x, w))dxF{x, w) 


{d 2 Y)ix, F{x,w)) {d 2 ip)ix,F{x,w)) 

{di'ip){x, F{x, w)){d 2 ip){x, F{x, w))dxF{x, w) 
{{d2'4’){x,F{x,w))f 


We denote each of the three terms on the right hand side above Di , D 2 , , respectively, and analyze them 

separately. For the first term, we have: 

Di{x,w) _ F{x,wY+^ J^^^^^^{dp){x,r)g{r)r-‘^-^^e{r)dr 


w 


u{x, F{x, w))g{F{x, w))4>e{F{x, w)) 


which can clearly be proved to be bounded using analogous steps to those used for dxF{x, w) after Eq. (A.3). 
For the second term, 

D 2 ix,w) _ ^{dii'){x,F{x,w)) f dxF{x,w) 


w u{x,F{x,w)) \ w 

which is also clearly bounded since \dxF{x,w)/w\ is bounded. For the third term, 

D3{x,w) _ {dxF{x,w)f dr{iy{x,r)g{r)^Yr)r-'^-^)\^^p^^^ 


w) 


w w i/{x, F{x,w))g{F{x,w))(j)s{F{x,w)) 

The above derivatives with respect to r generate the following simplified terms: 

{dxF{x, w))^ F(x, w)(d 2 k')(x, F{x, w)) {dxF{x, w))'^ F{x, w)g'{F{x, w)) 


w i'{x,F{x,w)) 

{dxF{x, w))'^ (t)'e{F{x, w)) 


W (f>e{F{x,w))' 


— {oi + 1) 


’ wF{x,w) 
dxF{x, w) 


g{F{x,w)) 


w 


All the above terms are clearly bounded in light of the Conditions (ST]|h) and (S4) and the facts that 
\dxF{x,w)/w\ and \dxF{x,w)/F{x,w)\ are bounded as already proved above. 


Finally, to show (4.8-iii), note that, there exists an e > 0 small enough such that for any \w\ < e, 

|1 + dx6{x,w)\ = |1 + duix, F{x,w)) + d 2 'y{x, F{x,'w))dxF{x,w)\ > g, 

for some r/ > 0, due to the Condition ( |S3| ) in Section and the fact that \dxF{x,w)\ < k\w\, for some 
constant k. ■ 


A. 2 Lemma 15.11 

Proof: We want to obtain the second order expansion in t of the function 
(A.4) H{t,z,q):=E u{z, J) F [Yt{e,(!},v) Y 
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where Y(s, 0, x) is the solution of (4.2), which has the same distribution law as X^(e,0,x). Let us start by 
representing (A.4) in terms of the density r(-;z) of z + 'y{z, J), the inverse of the mapping r —>■ j{z,r) := 
z + 'y{z, r), and the inverse of the diffeomorphism : r/ —>■ 0 , 77 ): 


(A.5) H{t,z,q):=E[ T{r-,z)iz{z,^ {z,r))l{Yt(e,ili,r);^q}dr ] = 




T{r;z)v{z,^ ^{z,r))dr 


Let Yt{q) := (q) and recall from the proof of Lemma C.l in [TO] that Yt{q) is a solution to the SDE 

pij pij ^ 

yt{r])=V- be{Yy{r]))dv + a'{Y^{r]))cr{Yy{r]))dv + j a{Yy{r]))dw'l + ^ S(Yy_(?]}, AZ*^}, 


/o 


0<v<t 


where S(u,() ■= d(u, —(} — u with d{u,() being the inverse of the mapping z ^ u := z + 6 {zX)- Above, 
{^7'lo^tsgr and are the time-reversal processes of W* and the Levy process Z^ := rp,* {ds, dr) 

(see, e.g., [221 Section VI.4] for information about time-reversibility). In particular, is a Wiener process, 
while Z*, is an independent Levy process with the same law of —Z* (cf. Theorems VI.20 & VI.21 in [22]). 
Next, we write (A.5) in the form: 

H{t,z,q)=¥.\HJYt{q 


where Hz{q) := r(r; z)i'{z,j ^{z,r))dr, and aim at applying the Dynkin’s formula (3.8) to deduce: 


(A. 6 ) H{t, z,q)=E { V(g)) = {q) + t{L,Hz) [q) + (1 - a)E [(L^i?.) ( Y^tiq 


da, 


where is the infinitesimal generator of {Yt{il)}t^o- We now proceed to justify (A.6) and the desired 

3 that the following t 

iu)LjGC!{R), for feC^( 


boundedness conditions (5.3). To this end, it is easy to see that the following two conditions are sufficient: 

(i) Hz G 


(A.7) 


The condition (A.7-i) follows directly from Condition (SI) and Lemma 3.1 For the other condition in (A.7), 
let us note that (see [TO] Eq. (3.5)]) 

(A. 8 ) 




Lef{y) ■■= be{y)f'{y) + -^f”{y) +T^f{y), where 


(A.9) 


^ef{y) = 


f'iy + b{y, r)/3)(l - p)d/3 ( 5{y, r)) hei-r)dr. 


and b^{y) := —bs{y) + a'{y)a{y). Clearly, the first two terms in L^f belong to (^^(M) when / G C^, provided 
that a'^,be G which holds true in light of Conditions (SI), (S2 and ( Sifi ). The boundedness of 2^/ 
would hold true provided that / G and 


(A.IO) 


sup|5(ri,r)| < fe|r|, r G (—e,e). 


for some constant k. Since 5{u,0) = 0, it suffices to show that sup^gjj \dr5{u,r)\ < 00 , where := 
{—£, 0) U (0, e). Note that 


drd{u,r) = -{d 26 ){u,-r) = 


{d 2 S){ 6 {u,-r),-r) 

1 -h (9i(5)((5(tt, -r), -r)'‘ 


r G De, 
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which is bounded on M x by Proposition 4.1 The formal differentiation of X^f{y) yields 

"1 . ^ ^ 2 


(A.ll) 


dyXsfiv) = 


'y^sj\yj-i I + hy^r)f3){l -/3) (l +{did){y,r)l3'^ d/3 (^6{y,rfj he{-r)dr 

f"{y + Hy, ^)/^)(l - /3)dl3S{y, r){did){y, r)hs{-r)dr, 

0 


+ 2 


and it is clear that for the derivative to be well-defined and bounded it suffices that / G CKM), (A. 10), and 
(A.12) 

Note that 


sup |(5i5)(tt, r)| < fc|r|, rG(—e,e), 

U 

{di6){6{u,-r),-r) 


idi5){u, r) = {di5){u, -r) - 1 = 


1 {diS){5{u,-r),-r)' 


r G D„ 


which can be bounded by k\r\ on M x by Proposition 4.1 


Formally differentiating (A.ll), generates the following terms: 

(A.13) dyief{y)= I ^ f^'^\y + 6{y,r)/3){l-P)(l + {di6){y,r)P^ dp (^6{y,r)^ he(-r)dr 

f^^\y + 5{y,r)P){l - P){d‘f6){y,r)Pdp (d{y,r)j he{-r)dr 


+ 

+ 


^ [ [ f^^Hy + Hy,r)P){l - P) (l + {di6){y,r)P^ dp6{y,r){did){y,r)he{-r)dr 


+ 2 
+ 2 


f{y + S{y, r)P){l - p)dp ( {di5){y, r)) he{-r)dr, 


f{y + S{y, r)P){l - P)dp5{y, r){df6){y, r)he{-r)dr. 


The previous expression shows that in order for the second derivative to be well-defined and bounded it 
suffices that both conditions (A.10) and ( A.12| ) are satisfied as well as 

(A.14) sup |((9i^)(u, r)| </i:|r|, rG(—e,e). 

U 

The latter condition again follows from Proposition |4.1| since, for r G D^, 

{df6)i6{u,-r),-r) 


{di6){u,r) = - 


(l -h {di5){6{u,-r),-r)y 


Some further algebra shows that the term s (g) and {L^Hz) (q) in (A.6) coincide with expressions of 
Ho{z]q) and Hi[z,q) given in (5.2). From (A.7) and the fact that the derivatives Hy\q), f = 0,... ,4, are 
bounded both on z and q, it is clear that sup^ g \ {L^Hz) (g) | < oo and sup^^g \ {L‘^Hz) {q) \ < oo, which in 
turn imply the conditions in Eq. (5.3). ■ 
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A.3 Lemma 15.21 


Proof: Throughout the proof, we write := li(e,0,z). Let us start by noting that the SDE (4.2) 
satisfies the Hypotheses 5-9 in [3], with E = (—e,e), G{dr) = ^^{r)h{r)dr, and rj{r) = r, since obviously 
/fg \r]{r)\Ph{r)dr < oo, for any p ^ 2, and, by Theorem 


4.1 


both 


Hz,r) 

r 


and 


dz6{z,r) 

r 


are bounded. Then, Theorem 5-10 in [H] guarantees the existence of a unique solution of (4.2). In turn the 


latter shows the existence of a unique weak solution for the SDE (3.3). Using an interlacing technique similar 


to that used in Theorem 6.2.9 of pQ, one can proceed to show the existence of a unique weak solution to 


(3.1) , which in turn implies the existence of a unique weak solution for the SDE (2.1). 

To conclude that x —)• Yf is a diffeomorphism, let us frame the family of solutions of the SDE 

(4.2) in the form (5-22) of [3], indexed by the initial state z in a bounded neighborhood U of x G M, with 
= z, = be, = a and = 5. First, note that the SDE (4.2) satisfies the conditions (i), (ii), 


and (iv) of [3l Hypothesis 5-23] since the coefficients be{y), (7{y) and 5{y,r)/r are twice differentiable in y 
and their respective partial derivatives with respect to y are bounded in light of Assumptions 
The assumption (hi) of [3l Hypothesis 5-23], and (i 


Proposition 4.1 


(S3) and 


satisfied since = be, B 
5.24 in [3j, the mapping : x i—>■ T)‘ 

SDE obtained by formal differentiation of (4.2): 


ii) in [31 Theorem 5-24] are trivially 
= a,C^ = 5 are deterministic functions independent of z. Therefore, by Theorem 
is differentiable and its derivative, dxYf is the unique solution of the 


I* s s 

(A. 15) dxY^^ = 1+ y, (Y:) dxY^dv + / u' {Y^) dxY^dW: + di6 (Y^-, r) d^Y^. jf {dv, dr), 

JQ Jo Jo J 

Note that the coefficients of (|4.2|) are deterministic functions instead of stochastic processes, so the terms 


dA,dB and dC in (5-25) of [3] are absent in (A.15) above. In particular, dxYf is given by 
dxYf = exp (^t-^ cf' {Y^)^ Jo / d*idv, dr)'^ , 


where U := /q (T^^) dx-|- jQa'{Y^)dW* + Jq J di6 (Y^_,r) j2*{dv,dr). Due to (4.8), it is clear that, a.s., 
dxYf 7 ^ 0, for all t. Hence, the mapping admits a differentiable inverse by the implicit function theorem, 
and, finally, is a diffeomorphism on M. ■ 


A.4 Theorem 16.II 


Proof: Throughout the proof, we only assume that v{x, r) ^ h{r) and the function v{x, r) 
admits an extension on M x M that is in x, which is weaker than the technical assumption 


El 


:= u{x,r)/h{r) 
(SI) in Section 


The case Nf = 0. 

Recalling that we denote the law of a process V (respectively, the conditional law of V given an event H) by 
E(V) (resp. C{V\B)), we have 
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Consequently, the inequality (3.4) implies that 


Xl {£,x)^x + y Nf = 0 


^ Ct 


N 


for any y > 0, where N > 0 can be made arbitrarily large by taking e > 0 small enough. 
The case Nf = 1. 


Conditioning on the time of the jump, we have 


Xf {e,x) x + y 


N! = l 


t 


Xt{e,{s},x) ^ x + y 


ds. 


Let J and U denote, respectively, a random variable with density function he(r) = hc{r)/Xe and an indepen¬ 
dent random variable with standard uniform distribution on (0,1). Then, conditioning on Xg -, the integrand 
above can be written as follows: 


^ x + y 


= E 

= E 


Xf_, (e, %,v)^x + y 


. V=X^. (e,0,a;)+7(X^_ (£,0,a;), j)6»(X^_ (e,0,x), J,(7) 


'■C/>p(X^_(£,0,a:),j)^ 


(e, %,v)'^x + y 


- V=X {£$,x) 


-hE 
= E 
+ E 


'-U<p{X^-{efi,x),jy 


X^_g (e, x + y 


- (£,0,x)+ 7(T^_ (£,0,3;),j)_ 


1 - 1? 


(^X^-(e,0,x), j)) E Xf_, (e,0,u) '^x + y 


- v=X 


V (x^-(e,0,x), J) E ^ x + y 


- 'v=X^-{£fi,x)+'y(x:^-(e,(/),x),j) 


where {Xs(e, 0, a:)}o^s^t is an independent copy of {Xf (e, 0, x)}o^s^i. We denote the last two terms above 
by Ti{x,y) and T 2 {x,y), respectively. By the fact that u ^ h and from the Markov’s property. 


Ti(t) ^E 


Xf_^ (e, 9,z)^ x + y^ 


Z=X^— {€,$,x) 


= E[Xf(e,0,x) '^x + y], 


which can be made 0{t^) for any y > 0 and an arbitrary large > 0 in light of (3.4). On the other hand, 

T 2 {x, y) = e\h (t - s, Xg- (e, 0 , x),x + y] 


where 


H{t, z,q) =E 


u{z,J)¥ Xf{£,%,v)^q 


- v=z+^{z,J) 


Using Theorem 5.1, T 2 {x,y) can be written as 

(A.16) r 2 (x,y)=E Ho (^Xg-{£,^,x);x + y'^ + {t - s)Hi (^Xg-{£,II},x);x + y^ +0{f), 

ij). By writing Ho{z-,q) as A~^ f-^i(z,q-z) i^{z,r)he{r)dr and 


where Ho{z]q) and Hi{z]q) are given in (5.2 
recalling from Lemma 3.1 that 7 “^ is and t 


le regularity on v imposed by (SI), it follows that Ho{z\ q) is 


in z. Then we can apply the second order Dynkin’s formula (3.8) to E[Ho{Xs{£,^,x)]x -by)] to get 


(A.17) 


E 


Ho [Xg{£,^,x)-x + y] = Ho,oix-,y) -b siLo,i(®;y) + C>(s ), 
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where 


1 

(A.18) Hofl{x-,y) := Ho{x-,x+ y) = — Ve{x,r)dr, 

J'y~^(x^y) 


(T^(x) 


(A.19) y) ■= LeHo{x] x + y) = bsix)diHo{x] x + y) + ^£^^dlHo{x] x + y) + Ho^i{x\ y), 

with 

HQ,i{x]y) = j [Ho{x + 'y{x,r)]x + y) - HQ{x]x + y) - diHo{x]x + y)'y{x,r)]i'eix,r)dr. 

From the dehnition of Hq{z] q) := A“^ f^i(z q-z) r)dr, one can readily check that 

roo 

(A.20) XediHo{x] X + y) = (x, 7 "^(x,y)) [9i7"^(x,y) - d 2 '^~^{x,y)\ + / diUe{x,r)dr, 

J'y-^{x,y) 

(A.21) \edlHQ{x]x + y) =-d2Ve{x,-^~^{x,y)) {di'^~^{x,y) - d2i~^{x,y)) +/ dlve{x,r)dr 

J 'y~^ix,v) 


-Ue{x,'y ^{x,y)) [91^7 ^{x,y) - 28182 ^ \x,y) + 8 ^j ^{x,y)] 


(A.22) 


XeHo,i{x-,y) = 



7~^(x+7(x,r),^—7(x,r)) 


i^e{x + 'y{x,r),ri)dri 


/ Ve{x, ri)dri - Xe 8 iHo{x; x + y)j{x, r) j Ve{x, r)dr. 

J'v-''-(x,y) J 


Next, in order to apply the hrst order Dynkin’s formula (3.7) to E[i7i(Xs(e, 0, x); x + y)], we need to first 
show that the mapping z 1 -^ Hi{z] x + y) belongs to C^. Since Hi{z-, q) = D{z-, q) + I{z] q) with D[-; q) G 
by the regularity of v{z,r) in z and by Lemma 3.1-1, we only need to show that I{-',q) G C^. To this end, 
let us write 


(A.23) 

where 


I 82 C{z,Pr;q){l - I3)d(3h{r)(j)s{r)r dr, 


C{z,r]q):= / ^(z,r))) 82 j ^{z,y)v{y,r)dy - ^{z,q)) 82 ^ ^{z,q)-f{q,r)i^{q,r) 

J'y{q,r) 

The representation ( A.23[ ) follows by the Taylor’s theorem and the fact that C{z,0-,q) = 82 C{z, 0 -,q) = 0. 
Next, observe that C{z,r;q) is in {z,r) since all the involved functions are C^. Therefore, standard 
applications of the dominated convergence theorem implies that I{z;q) is in 2 ; and, furthermore, 

dll{z-,q) = j 8 {{ 8 lC{z,ldr-q)){l-I3)dl3h{r)^e{r)r‘^dr, i = l,2. 


Then, we apply the first order Dynkin’s formula (3.7) to K[Hi{Xs{^, 0, x)] x + y)]: 

(A.24) E iLi (^A^-(e,0,x);x + y^ = iL 7 o(x; y) + 0(s), 

where, using the following relationships obtained by implicit function theorem 

7 “^(x,x + y) = 7 “^(x,y), 82 T^ix,x + y) = 82 'y~^{x,y), 5|7“^(x,x + y) = (9|7“^(x, y), 
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we have 


y) := Hi[x\ x + y) = Difi{x] y) + hflix] y), with 
Di,oix-,y) := Y(^[beix + y) - v'{x + y)]us {x,j~^{x,y)) d 27 ~^{x,y) 


^^25) -v{x + y)[d2i^e{x,'y ^{x,y))[d2'y (x,7 ^{x,y))dl-f ^(x,?/)]). 

hfi{x-,y) ■= Y j / i'6{x,T^ix,y)) d 2 T^{x,y)ueiv,r)dy 

'^e J _ J'y(x+v,r) 


^x+y 

. J^{x+y,r) 

- {x,'y~^{x,y)) d2'y~^{x, y)-/{x + y,r)ue{x + y,r) 

Finally, substitute (A. 17) and ( |A.24[ ) into (A. 16), we have 


dr. 


(A.26) 


(e> {■s}, x)'^ x + y = Hoflix] y) + sHq^i{x; y) + {t - s)Hifl{x] y) + O(t^), 


where 77 o, 0 )-f^o,i and are given in ( A.18[), (A.19) and (A.25), respectively. Thus, 


(A.27) 


Xf (e, x) ^ X + y 


Nt = l 


= Ho,o{x-,y) + -[Ho^i{x-y) + i 7 i,o(x;y)] + 0 {f). 


The case Nf = 2 


Conditioning on the times of the jumps, we have 

2 rt rt 


(A.28) 


A'i(e,x) ^ x + y 


N! = 2 


t 2 


0 J Si 


Xf(e,{si,S2},x) ^ x + y 


ds2dsi. 


Denote the probability inside the integral above by ^^(si, S 2 , x, y). Below, let Ji and Ui {i = 1,2) be 
independent copies of J and U, respectively. Denoting an independent copy of {Xf (e, {si}, x)}o^s^i by 
{Xs(e, {si}, x)}osgssgi and conditioning on X^-, we have 


(A.29) 

where 


At{si, S2,x,y) = E 
= E 


Xf_ (e, 0 , x) ^ X + y 


'^S2,t (X^2(e>{'Si},a; 


- v=X _ (£,{si},x)+7( X _ (£,{si},a:),J )6»( X _ (e,{si},x),J2,!72 


+ E 


'^S2,t (X^ 2 (e>{si}x 


^S2,t(a^i) := IE[1 - i^(xi, J2)]E Xf.^^ (e, 0, xi) ^ x + y = E[(l - i/(xi, J2)) (e$,xi)'^x+y]^ 


Ts2,t(xi) := E 


V (xi, J2) E Xf_ (e, 0 , x) ^ X + y 


v=xi+'i(xi,J 2 ) 


We denote the last two terms on the right-hand side of (A.29) by T^{x,y) and T 4 (x,y), respectively. Condi¬ 
tioning on X -, we can write T 3 as 


T3 = E 
= E 


E 




2 ,i (Xs2-si(e>0X2)) 


X 2 =Xs^ (£$,x)+'i(X,^ (e$,x),Ji)e(X,^ (£$,x),Ji,Ui) 

l-u[Xs^ (e, 0 , x), Ji) ) E Ts2,i (x52-^i (^) 0 : X2 


(A.30) 


-F E 


= : T37 -I- r 3 , 2 , 


V (Xs^(e, 0 ,x), Ji) E Ts2,t (X 52 _si(e, 0 ,X 2 


a: 2 =Xsj {e$,x) 

X 2 =Xs^ (£,0,3;)+7(Asi (£,0,a:), Ji)_ 
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where {Xs(e, 0, x)}o<gs<gi and {^^(e, 0, x)}o<gs^t are independent copies of {X®(e, 0, x)}o<gs<gt. Using that 
0 ^ ^ 1 and the Markov’s property, it is clear that 


173,11 

= E 
^ E 


E 


(A.31) 


^S2,t —Si (S; 0) 2:2)^ 

'l/(x,,(e,0,x))| 

^ x + y 

Xf {e, 0, x) ^ X + yj = 0{t), 


X 2 =Xs^ {e,(D,x) 


xi=Xs„{e,<D,x) 


as t —)• 0, by taking e > 0 small enough. To deal with the second term in (A.30), let us define 
'^S2,4(3^1) •— E ^1 — i/ '^{Xl_^^{efli,x\):iiX+y\ J 

and note that 

$,2,i(xi) - T,2 ,i(xi) = E (^P(xi,J2)-i^(xf_,2(e,0,xi),J2)) 

^ sup |9ii/(x, r)|E xf_^^{£,%,xi) - xi 
x^r L 

Using the facts that Xf_g^ (e,0,xi) — xi = Xf_g^ (e, 0,0) and sup^,^ |9iz^(x,r)| < 00, the last term above 
converges to 0 as S2 < t —)• 0. Therefore, it suffices to study the asymptotic behavior of the expression 


73,2 : = E 
= E 


(^Xsi(e,0,x), Ji) E T^2,t (X^2 -si(£, 0,3^2 

P fXsi(e,0,x), Ji 


X 2 =Xs^ (£,0,x)+'y(Xsj^ (e,0,x),Ji)_ 


X E 


1 - 


^ Si (^) 0) ^ 2 ) ’ '^2 j j (£,0,X2)^x+y} 


X2=Xb^ (£,0,x)+7(Xsj {£,0,x),Ji) 


Using again that sup,j,,, \y{x, r)| < 00, sup^, \dii'{x, r)| < 00 and E [|X^ (e, 0, x) — x|] = E [|Xf (e, 0, 
0, uniformly in x as s —)• 0, 


(A.32) 

where 


73,2 = E 


i7 (Xsj(e,0,x),t - si,x + y) +o(l). 


H{z, t, q) := E 


X 2 =z+'y{z,Ji) 


'{z,Ji) |(1 - i/(x 2,72))E xf (e,0,X2) ^ y | 

Using a similar argument as in the proof of Lemma |5.1[ 

(A.33) H{z, t,q) =E [v (z, Ji) {1 - i? {z + 7(2, Ji), J2)) lz+'r{z,Ji)^q] + tTZliz; y), 


where sup^ ^ \TZj{z‘, y)| < 00 for t small enough. Furthermore, the first term on the right-hand side of (A.33) 
can be expressed as 

r TOO ^ 

u{z,j~^{z,r)){l - i>{r, r 2 ))T{r; z)drh£{r 2 )dr 2 , 
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which is C? in z. The previous fact together with Dynkin’s formula as well as (A.32)-(A.33) imply that 


^3,2 = IE 


i^{z, 7 (z, r))(l - z^(r, r 2 ))r(r; z)drh^{r 2 )dr 2 


x+y 


z=Xsi (e$,x)_ 


+ o{l) 


v{x, 7 ^(x, r))(l - z/(r, r 2 ))r(r; x)drhe{r 2 )dr 2 + o(l). 


x+y 


Using the change of variable ri = 7 ^(x,r) and the representation 
(A.34) f(y;z) = dy [ lz+^(z,r)<yKir)dr = ^ d2T\z,y), 


A. 


together with ( |A.30[ )-( |A.3l[ ), we have 

T 3 = ^ / / i/(x,ri)(l - z^( 7 (x,ri),r 2 ))/i£(ri)(iri/i£(r 2 )(ir 2 + 0 ( 1 ) 

rOQ poo 


(A.35) 


i^eix,ri)dri - ^ I iys{x,ri) i^eilix,ri),r2)dr2dri + o{l). 

'^£ J'y~^(x,y) '^e J'y~^(x,y) J 


We next consider the second term, Ti{x,y), in (A.29). By Lemma 5.1 


T 4 (x,y)=E Ho[Xs2i£,{si},x)]x + y) +0{t) 


Conditioning on X -, 

®i 


T4{x,y) = E 

E 

E 

Ho (£) 0 w); + 2 /) 






v=z+"f{z,Ji)d{z,Ji,Ui)_ 

z=Xs^{s,0,x)_ 


+ o{t) 


= E 


E 


(1 - u{z, Ji))E Ho (J(s2-si (e, ^,z)',x + y'^ 


z=Xsi {efi,x)_ 


(A.36) 

(A.37) 


+ E 

= E 

+ E 


E 

1 ^( 2 , Ji)E 


- 


- v=z+'f{z,Ji) 

E [(1 - i?{z, Ji))] E Ho (Xs 2 -si (e, ^,z);x + y 


Z = Xs-^^ (£,0,x)_ 


+ 0{t) 


z=Xsi (eM,x)_ 


iy(z,j (z,r))E Ho I Xs^-sAsJh}; x + y) r(r;z)dr 


z=Xs^ (e,0,x) 


+ 0(t) 


Ho (Xs 2 -si(£,<^,z);q) -Ho(z;q) ^ sup |(9^iLo(2:w)|E Xs^-s^{e.,^, z) - z 

^ ^ z,q L 


Let us denote the expressions in (A.36) and (A.37) by r 4 ^i and T 4 ^ 2 , respectively. Next, since Ho{z\q) is C^, 
(A.38) E 

L \ / J z,q 

'' D ^ 

Using the facts that sup^^^ \dzHo{z;q)\ < 00 and (e, 0 , 2 ;) — z = Xs 2 -si (£,0,0), the last term above 

converges to 0 uniformly in 2 ; as S 2 — si < t —0. Therefore, recalling the definition of Ho{z;q) = 
f^Hz,q-z) MzH)drlXe, we have 

T4,i = E 


E[(l -J 2 ( 2 , (A^i(£,0,x);x + yj + o(l) 


(A.39) 




poo p 

/ v^{z,r)dr / {1 — u{z,r)) h^{r)dr 


J 'y~^{z,x-\-y—z) J 

z=Xs^ {e,0,x)_ 


+ 0 ( 1 ). 
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Since the expression inside the expectation above 


in z, using again that E 


Xs-^ (e, 0 ,x) - X 


las partial derivative with respect to z, uniformly bounded 


= E 


(e,0,O) 


0 uniformly in x as si < f —)• 0, 


1 


^4,1 = V 


Vs{x-,r)dr / (1 — i/(x, r))/ie(r)dr + o(l) 


(A.40) 


J'y 1 {x,y) 

I V£{x,r)dr - ^ I Ueix,r)dr j iy£{x,r)dr + o{l). 


We next consider the term Z 4 ^ 2 (i) in (A.37). Using (A.38) again and an argument similar to that from 

+ 0(1) 


(|A.39|) to (|A.40p, we deduce as follows. 

v{z, r))Ho{r-x + y)r(r; z)dr 

/ poo 

P(x,7"^(x,r)) J ^ 


Ta2 — IE 


1 


z=Xs^ (£,0,3:)_ 

Veix, r 2 )dr 2 T{r\x)dr + o(l) 


/7 ^ {r,x-\-y—r) 

(A.41) ^ +2 / / U£{x,r 2 )dr 2 dri + o{l), 

where, in the last equ ality, we us e the change of variable ri = 7 “^(x,r) and the representation (A.34) of 
T{y;z). Summing up (A.40) and (A.41), we have 


(A.42) 


roo 2 r'^ r 

Ty^{x,y) = — / Ve{x,r)dr - ^ / U£{x,r)dr / U£{x,r)dr 

'^e J'y~^(x,y) '^e J'y~^(x,'ij) J 


e Jz~^{x,y) 


+ 4/„.(.,r.)/ 


7-'-(a;+7(x,ri),j/-7(a;,ri)) 


V£{x,r 2 )dr 2 dri + o(l). 


Therefore, summing up T^{x, y) and T 4 (x, y) from (A.35) and (A.42), P [Xy{e, {si, S 2 },x) x + y] = H 2 fl{x] y)+ 
o(l) with 




2 1 f 

fi{x]y) = -^l i^£{x,r)dr-—^ iy£{x,ri) i^£{^{x,ri),r2)dr2dri 

J'y~^{x,y) J'y~^(x,v) J 

coo 


e 47-1(3;,i/) 


(A.43) 


poo p 

/ V£{x,r)dr / i>£{x,r)dr 

47 - 1 ( 3 ;,I/) 4 


\ 2 

47-1(3;,3/) 


1 

+ Af 


/ poo 

Ue{x,ri) J 


U£{x,r2)dr2dri, 


and, due to (A.28), 
(A.44) 


■y-''-(x+-f{x,ri),y--y{x,ri)) 


Af(e,x) ^ x + y Nf = 2 = H 2 fl{x]y) + o(l). 


Substituting (A.27) and (A.44) into (6.1) and using the fact that e +* = 1 — —h 0{t^) as t J, 0, 


we obtain 


where 


¥[Xt x + y] = tPi(x,y) + —P 2 {x,y) + 0(t^). 


poo 

(x,y) = XeHofiix-,y) = / i^e{x,r)dr, 

J-y-^{x,y) 

P 2 {x, y) = -2\lHofl{x\y) + Ae[i7o,i(a;; v) + Hifl{x; y)] + X‘^H 2 ,o{x; y). 


P 
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Some further simplifications lead to the expressions of Pi{x,y) and P 2 {x,y) stated in the statement of the 
theorem. ■ 


A.5 Corollary 7.1| 

Proof: Let us hrst note that, due to the Condition (S3 i]), 

| 7 (x,r)| ^ cr, 


for all X and r, where c is defined as in Condition (S5). Next, define g^{r) := g{r)e^'' and note that, in view 


of the Condition (S5), h^{r) := g"^{r)\r\ “ ^ is a valid Levy density. Furthermore, < 7 ^ clearly satishes the 


other requirements of Condition (S4), while 


_e v[x,r), 


can be readily seen to meet all the requirements of Condition (SI). The result is then a consequence of 
Theorem 16.11 ■ 


A.6 Lemma 18.11 


Proof. From the infinitesimal generators (8.2)-(^8.3L we can identify the semimartingale characteristics 


{B, C, -i?) and ,'d^) of {Xt}t^o and respectively, relative to the truncation function 

(A.45) Bt = f b{Xs)ds, Ct = f a‘^{Xs)ds, 'd{dt,dr) = K {X^-,dr) dt, 

Jo Jo 

pi pi 

(A.46) ^ y ^ ^ y ^ 

where K {x, A) := f 1 a (7 (x, r)) r' (x, r) dr, b^{x). We also dehne bs{x) and (x) as 

(A.47) beix) := b{x) + j ^^(x, = b{x) + j 7 (x, r)l||^(,,_^)|>i}i/(x, 

(A.48) al{x) := al{x) + j K^ix, dr)r‘^ = a‘^{x) + a'^{x) + j 7 ^(x, r)i/(x, r)l{|^( 3 ,_^)l>£}dr. 

Clearly, for e G (0,1), be = b{x) + f K{x,dr)rli^\.e\>i} and d1{x) = ct^(x) + / K{x,dr)r‘^, in light of the 
definition of dg. Also, 

I Keix,dr)gir) = J < 7 ( 7 (x, r))Kx, r)l||,(,,,)|>,|dr J gi^ix,r))uix,r)dr = J K{x,dr)g{r), 
for any bounded continuous function g that vanishes in a neighborhood of the origin. Therefore, by m 

~ 'jy i-r-w 

Theorem IX.4.15], {Xf}t^Q — {Xt}t^o provided that the uniqueness hypothesis [151 IX.4.3] holds true. By 
na Theorem III.2.26], the uniqueness requirement stated in |15l IX.4.3] is equivalent to the weak uniqueness 
of the solution for the SDE defining X. This fact was established in Lemma |5.2[ Therefore, we conclude the 
convergence result claimed in the lemma by |151 Theorem IX.4.15]. ■ 
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